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Abstract. Wo develop an effective field theory for lattice models, in which the 
only non-vanishing diagrams exactly reproduce the topology of the lattice. The 
Bethe-Peierls approximation appears naturally as the saddle point approximation. 
The corrections to the saddle-point result can be obtained systematically. We 
calculate the lowest loop corrections for magnetisation and correlation function. 
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There are two ways for doing mean-filcd-likc approximations in lattice models. 
The first one consists in replacing interaction constants decaying with distance by 
effective distance-independent interactions. The best known of these approximation 
schemes is the Bragg- Williams approximation P0|. The second path to a mean-field 
result neglects the correlations induced by the presence of closed loops on the lattice. 
This is the Bethe-Peierls (BP) approximation 0. The former approach is exact in 
fully connected models, while latter is exact on the Bethe lattice jH). 

The list applications of the BP approximation in statistical physics is very long. 
As an example, let us mention only the application to spin glasses [7] and 
combinatorial optimisation jS], which attracted much attention recently. Within this 
context, the BP approximation essentially coincides with the cavity method. 

The systematic field-theoretical treatment of the corrections to the first method is 
well developed and widely used. On the other hand, to our knowledge, the corrections 
to BP results were not yet formulated in a similar manner. The actually available 
approaches rely mostly on cluster variational methods [S] or expansions in inverse 
dimension |10j . where loops of limited size arc taken into account. Recently, a very 
promising path was paved by Montanari and Rizzo who were able to consider 
loops of arbitrary size, correcting the BP result for fsing model on hypercubic lattice 
and also improving the cavity equations for a spin glass. 

An important step the direction of field-theory formulation of the BP 
approximation was taken by Efetov |I2| in the context of supersymmetric formulation 
of the Anderson localisation problem. The method of Efetov consists in writing an 
effective Lagrangian, with the property, that the sequence of Feynman diagrams of 
the corresponding field theory is equal to the sum of partition function of the original 
model, plus partition function of the same model wrapped several times around. It is 
expected that in thermodynamic limit the wrapping is unessential. 

Then, the saddle-point equations for the effective field theory correspond exactly 
to the BP approximation. In this way, Efetov was able to locate the localisation 
transition and calculate the correlation function to the zero-loop order. It is worth 
noting, however, that the zero-loop order on the correlation function corresponds 
already to the lowest correction to the BP result, in which the connected correlation 
function is exactly zero. 

The purpose of our present work is to show, how the Efetov's method can be 
applied to general lattice model, especially magnetic models like Ising or Ginsburg- 
Landau models. We show, how the well-known BP results are generated within 
the formalism and calculate the lowest-order corrections, which are one-loop for the 
magnetisation and zero-loop for the correlation function. 

We shall consider a system defined on a lattice represented by an oriented graph 
{M, £}, where M is the set of nodes, while £ C Af x AT is the set of oriented 
edges. For each node r G Af define T r = {a G £\3r' G Af : a = (r',r)} and 
O r = {a G £\ 3r' G Af : a = (r, r')} the sets of incoming and outgoing edges, 
respectively. For further convenience we denote a the initial and a the final node 
of the edge a G £ , i. e. a = (a, a). 

Each node r is described by the state variable <f> r G S. We shall treat both discrete 
and continuous 4> on equal footing, denoting by summation in discrete case and 
integration in continuous case. The state of the system is described by $ G . For 
the Efetov method to work it is essential to decompose the Hamiltonian of the system 



Loop expansion around the Bethe-Peierls approximation for lattice models 



3 



into node and edge parts, 

H= Hx + H £ 

H N = V A<Pr) (1) 

reA/ - 

r££ 

Therefore, all the terms in the partition function are products of node and edge factors, 



S> $ \a££ / VrGJV / 



(2) 



where j((f>,<f>') = e u ^'^ > and Cr(<^>) = e^^'. When appropriate we will use also 
matrix notation ( r ((j>,(f)') = ( r (<fi)5((j) — </>'). We will also use the inverse matrix 7 -1 , 
where obviously £ „ 7 (0, ft'h'Ht", = H<t> ~ <t>') 

In order to build the effective field theory, we introduce complex auxiliary fields 
for each edge a £ £ and each value of the physical field (j> £ S, x a (<f>) G C and its 
complex conjugate x a {4>)- 

The Efetov Lagrangian is similarly decomposed into edge and node part 

£ = £o + £i (3) 
where the edge part is simply 

£o = EE S-(07- X fa, 0'K(0') (4) 

a££ ct>ct>' 

while the node part is slightly more complicated 

a = - x; x ( n ( n ) ■ ( 5 ) 

The key quantity is the Efetov partition function 



x,x\e 



-C[x,x] 



Z= fd\x'i\e-*>M ={e ~ Cl)o (6) 
where we denoted the average with respect to Co as { )q. 

The diagrams contributing to © correspond to the lattice of the original system, 
wrapped m-times around, where m = 0,1,2,3,.... Indeed, the first contribution is 
obviously 1 and in all the remaining diagrams the lines are exactly the edges of the 
lattice and the graph vertices are the nodes. Each node and every edge must appear 
the same number (which is m) of times, but they may be connected by various different 
manners. That is how the wrappings come around. The linkcd-cluster theorem enables 
us to write 

Z = e z < (7) 

where Z c is the sum of all connected components coming from the wrapping. The first 
term is Z, second term contains various possibilities to wrap system twice without 
leaving two disconnected parts etc. 

The averages of physical quantities, like magnetisation and correlation function, 
arc calculated by making derivatives of Z with respect to Cr (</")• Using the notation 
d r {4>) = Cr(<P) dc^irj}) we ^ nc ^ that the magnetisation on site r is 
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and the connected correlation function C(r,r') can be obtained by derivative with 
respect to local magnetic field, giving 

Technically, the treatment is complicated by very high symmetry of the 
Lagrangian J2J), which cannot be easily parametrised. Let r\ a = e le = £ U(l) for each 
edge a £ £ . The field variables transform as 

X a {4>) -> r]aXa{<P) (10) 
Xa{<t>) -*■ VaX a (<P) (H) 

with one condition at each vertex (Jlaei 'faXlIaeo = wmcn can be 

interpreted as the requirement of current conservation . The number of free parameters 
of the symmetry group is \£\ — |7V|, i. e. proportional to the system size. 

As a zeroth approximation we can calculate the observable quantities using saddle- 
point method. This amounts to solving the coupled set of equations 

<9£[a;«,z«] _ dClx^^x^} _ Q ^ 12 ^ 



which means 



dx a (4>) dx a {4>) 

" (s ^)=z 7(^00 n 4 s V) n ^ s V) 



0' bela fceOa(6#a) 

si s) (0)= n 4*V) n 4 s) (0o- 

This set of equations has infinitely many equivalent solutions related by the symmetry 
group operations. In the following, we shall always select the real uniform solution 
Xa = %a = x\/a £ £ . 

To take a concrete example for the geometry, let us consider the c?-dimensional 
hypercubic lattice. We will use indices a, (3, ... £ {1,2, ...,d} to denote one of the 
Cartesian axes and indices a, 6,... £ { — 1,1} to indicate one if the two orientations 
along one of the d axes. The 2c? possible directions will be denoted by indices 
[i^v,... £ {l,...,2c?}. On hypercubic lattice it is convenient to change the notation 
for the field variables as x r 2a-i{4>) = x (r-e a .r)(4>) an d Xr2a(4>) = £(r.r+e„) (0)j where 
e a are the unit translation vectors along the Cartesian axes. 

The real uniform saddle point is the solution of |<S| coupled equations 

a#)=E^^W> 2d ~V)- (13) 



M = J TT^uHR (14) 



If we solve them, we can calculate the magnetisation 

E C(0)x 2d (0) 

while the connected correlation function is identically zero. These results coincide 
with the Bethc-Pcierls approximation, or, cquivalently, with the solution of the 
corresponding model on Bethe lattice with coordination number 2c? 
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Fluctuations around the saddle point are denoted as y ri i(4>) = x rtl {(p) — x(<j>) and 
the Lagrangian is 

C = C[X^ S \ X^] + \ J2 E E Mr^.M^ 0') Vr>utf) + - ■ ■ (15) 

rr' /ii-' 0^>' 

The Fourier transform of the bare propagator G = M _1 can be written explicitly 

G^(p|&^)= (16) 
= ^x 1 - d ^)C- 1/2 ^)B^,A)^ i/ (p|A) J B(0',A)C- 1/2 ^ / )^ 1 " d (^) 

A 

We must spend some time explaining the factors entering the expression l|16(l . 
The matrix function <?(p|w) = [w _1 mi(p) — mo] -1 with mo M „ = 1 — and 
iTT-iip)^ = ^a/3(l — S a g) depends only on the geometry of the lattice, but not on 
the nature of the spin variables 4>. Straightforward algebra gives 

S a p(l - S aS ) e ip,i u - 

(17) 



1-Lu 2 



2 uj- (c 1 ^ - w)(e- I P" -tu) 
~ » uUJ + l-2c(p)oj + {2d-l)ujZ 

where c(p) = J2 a =i cosp a . 

The numbers A form the set of eigenvalues of the (symmetric) matrix 

7(6 (/>') = .x d - 1 (^)C 1/2 (0) 7 (^ ^OC 172 ^')^" 1 ^') (18) 
and B((f>, A) are the corresponding normalised eigenvectors. Contrary to the matrix 
g(p\uj), the eigenvalues A and the eigenvectors bear the information on the nature 
of the variables </>, but conversely the dependence on the lattice structure is only 
weak. In fact, the lattice enters only through x(<p), which for transitionally invariant 
lattices depends only on the coordination number, as can be seen from the saddle-point 
equations l(T3|) . In such a way the Efetov formalism decouples to large extent the effects 
of the lattice structure from the effects due to specific choice of the Hamiltonian living 
on the lattice. 

It is easy to see that one of the eigenvalues of 7 is always equal to 1. Indeed, from 
the saddle-point equation multiplying both sides by x d_1 (0X 1//2 (</>) we conclude 
that the vector 

v 1 ^) = x d {4>)e /2 {(p) (19) 

is eigenvector of 7 and the corresponding eigenvalue is 1. This brings about a serious 
problem in the calculation of the diagrams, because the bare propagator g(p\uj) 
diverges for w — ► 1, as can be seen in Eq. (|17(l . The source of the divergence is 
the continuous degeneracy of the saddle point, as mentioned above. Usually such 
situation is solved by explicit integration over the degenerate saddle points, but 
here it is impractical, because the symmetry group cannot be easily parametrised. 
Fortunately enough, the physical observable quantities contain such combinations of 
the propagators that the divergences stemming from various terms cancel each other. 
We verified this property by explicit calculations of the magnetisation and correlation 
function up to one-loop level. 

In the diagrammatic expansions all summations over field variables <f> are reduced 
to sums over the set of eigenvalues A, so the lines correspond to the propagators g(p\X). 
The vertices are of two kinds. As a small filled circle we depict a vertex involving mere 
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summation over the field variable, i. c. • = J^, . . ., where the empty circle means 
summation of terms containing factor <j), i. e. o = </>.... Later, we shall see also a 
vertex denoted by empty square, denoting summation with factor <f> 2 . Generally, we 
denote the n-lcgged vertices as 

"-™D fe (Ai,A 2 ,...,A„)= (20) 



= vi 



lJljBfa, A,) 



1=1 



The analogous cases of empty circles and squares are distinguished by writing one or 
two overlines over D, respectively. 

Having all elements at hand, we proceed to computation of the one-loop correction 
to the magnetisation. We obtain pictorially 



M 





(21) 



Writing explicitly the value of the diagrams we can see, that the formal expansion 
parameter is ||vi|| -2 . We should note, however, that this parameter is not actually 
small; its value is at least equal 1. Collecting all the terms to the order ||vi||~ 2 we 
have finally 

M = ° 0+ nhp E [(D 2 (X,X)-D ) + 



1 

l --d 



E 



^ 1 (A') J D 3 (A',A,A)5..(0|A') g A {Q\X) 



(22) 



+ o(n« 1 ir 4 ) 

The combinations of the propagators occurring in the formula are 

2dA' 



3..(0|A')=Es^(0|A') 



1 - (2d - 1)A' 



and 



9A 



(0|A) ee £(1 - V)^(0|A) = 1 V d _V X2 
d d p 1 



1 + A 



* (2tt) c 



1 - (2d-l)A , 
We can see that the phase transition is located at A = A c 



(23) 



(24) 



(25) 



cosp Q )_ 

l/(2d- 1), as in the BP 
approximation. To calculate the correction to the transition temperature we would 
need to make partial resummation of the series, which is feasible, but the results will 
be presented in subsequent publication. Instead, we shall estimate the correction to 
the transition temperature by a simple argument later. 

As a special case, let us consider the Ising model. The matrix 7 has only two 
eigenvalues, one of which is always equal 1, as already stated. For now, denote A 
the other one. There is only single term in both summations of 1)221) and we have 
(superscript (0, 1) indicates that we include only zero and one-loop contributions) 



A/ (0 ' 1} = (cos2tf) 



1 



2IHP 



1 



1-^) 



(26) 
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( cos-d \ 

where the angle i9 is related by equation Vi = \\vi\\ ( g - n ^ ) to the eigenvector of 7 

corresponding to eigenvalue 1. This expression suggests that the applicability of the 
expansion relies on the smallncss of the effective expansion parameter </a/||ui || 2 . 

Far from the transition point we have A <C A c . Indeed, for small h we obtain in 
the high-temperature phase A ~ J, while in the low-temperature phase A ~ c - 4J ( d - 1 ) i 
so we can expand the result in powers of A. Keeping only the lowest terms, i. e. up 
to order 0(A 5 ), we get 

M ~ (cos 20) [l - t\ ~~T7"ir (1 + 2(d - 1) A) 4d(d - l)A 4 j . (27) 
L ll^l^ J 

In the saddle-point approximation the connected correlation function vanishes. 
However, the first correction is found already at zero-loop level. To be consistent with 
the calculation of the magnetisation, we proceed by differentiating the magnetisation 
with respect to local magnetic field. We find, diagrammatically (superscript 
indicating that we keep only zero-loop terms) 



C(»)( r -/) = i-Mf] 



^-[(. „- „) + (□. -oo)S(r-r')] 



(28) 



and therefore 

—2 



Do -Do" V D\(\) 



d d p e ipr (l + A) 



(29) 



-* W l -(2d- 1)A + ^ ££(1 - cosp Q ) ' 

Note that derivative with respect to magnetic field produces also on-site terms 
(proportional to S(r — r')) which contain summing over <f) a function multiplied once 

more by <f>. From here the new vertex n = Do = Y^<h 4" 2 B 2 ((I>, 1). 

For Ising model in high temperature phase at zero magnetic field Di(A) = — 1, 

Do = 0, Do = Do = 1 and we have 

C<°> (r) = f (30) 
J_ n (2n)d 1 - (2d - 1)A + ^ ^(1 - cos Pa ) 

We were able to compute also the one-loop contributions to the connected correlation 
function, verifying explicitly that the divergences stemming from the propagators 
g^ tL ,(p,uj) 1 ui — y 1 exactly cancel, leaving a finite result. The final formulae are rather 
lengthy, though, and will be reported elsewhere [T3] . 

Finally, let us consider the shift in the transition temperature for the Ising 
model. A consistent calculation would require explicit partial resummation of 
the loop expansion, as noted before. However, we can get an estimate of the 
correction by the following simple consideration. The magnetisation (|26(l has the 
form A/^ ' 1 ) = M^(h) [1 - m^(h)], where Af(°) is the BP (or saddle-point) value 
for the magnetisation and rrS^{h) is an even function of h, so lim/^o tt m Q 1 ) = 0- 
This follows from the fact that both A and ||wi|| 2 are even functions of h. 
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The expression (|26[1 can serve as a starting point for calculating first correction in 
the loop expansion for the inverse susceptibility x 1 — liniM^o g^j- So, up to order 
H^i || -2 we have 

(X- 1 )^ = (X^ 1 ) (0) [1 + ^ 1 )(0)] (31) 

and the shifted transition temperature is given by the condition m^^(0) = — 1. We 
get the following equation for A c = tanh J c 

A 1 , 2(d-l) Xc9a(0\X c ) (32) 

2d-l 2d- 1 2 ~r* (i _ \2)2^2 + 5A (0|A C ) ' 

This equation can be solved by iteration, the first step consisting in inserting the 
BP value A c 0) = (2d - l)" 1 to the right-hand side. This way we get the following 
estimates for the corrected transition temperature (the numerical integration was 
performed in Maple): J c = 0.2372651... for d = 3, J c = 0.1520003... for d = 4, 
and J c = 0.1144870 . . . for d = 5. These values are strikingly close to those calculated 
in Ref. |11| . Moreover, if we linearise the RHS of l)32|l in the supposedly small value 
of (^(OlAc), we get nearly identical expression for the correction as reported in |llj 
(equation (3.17)), the difference being only in a numerical factor which goes to 1 for 
d — ► oo. It suggests that our approach and that of Ref. JT] are intimately related. 
The details if that relation are not yet clear, though. 

To conclude, we developed a systematic procedure how to calculate corrections to 
Bethe-Peierls approximation, which can be used for any equilibrium statistical model 
on lattice. We showed it in the simplest case of the Ising model on hypercubic lattice, 
obtaining finally corrections to the transition temperature. In principle, our method 
should work on lattices of arbitrary geometry, but the practical calculation of the 
diagrams relies on translational symmetry An application to topologically disordered 
systems would be possible, but it would require some additional approximations. 
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